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Abstract 

Interwall interaction energies, as well as barriers to relative sliding of the walls along the nanotube 
axis, are first calculated for pairs of both armchair or both zigzag adjacent walls of carbon nanotubes with 
a wide range of radiuses. It is found that for the pairs with the radius of the outer wall greater than 5 nm 
both the interwall interaction energy and barriers to the relative sliding per one atom of the outer wall 
only slightly depends on the wall radius. A wide set of the measurable physical quantities determined by 
these barriers are estimated as a function of the wall radius: shear strengths and diffusion coefficients for 
relative sliding of the walls along the axis, as well as frequencies of relative axial oscillations of the walls. 
For nonreversible telescopic extension of the walls, maximum overlap of the walls for which threshold 
static friction forces are greater than capillary forces is estimated. Possibility of experimental verification 
of the calculated barriers by measurements of the estimated physical quantities is discussed. 

1 Introduction 

The possibility of relative motion of the walls [ll [2] in multi-walled carbon nanotubes (MWNTs) makes it 
promising to use nanotube walls as movable elements of nanoelectromechanical systems (NEMS) (see [3] 
for a review). At present, nanomotors in which the walls of a MWNT play the role of the shaft and the 
bush [H m El E] and memory cells operating on relative sliding of the walls along the nanotube axis [U [9] 
have been implemented. A number of NEMS have been proposed, which operate on the relative motion 
of carbon nanotube walls, including a gigahertz oscillator [10^ lllj . Brownian nanomotor [12], a bolt/nut 
pair \13\ HM [T5] and an accelerometer \16\ [T7] . Thus the study of relative motion and interaction of carbon 
nanotube walls is an actual problem. 

However neither the interaction between carbon nanotube walls nor the interaction between graphene 
layers have been investigated in details. The measurements of the energy of interaction between graphite 
layers give the wide scatter of the results (see [IHIIIHIEO] and references in [E]). Though experimental value 
of the interwall interaction energy is available (23 — 33 meV per atom [2] ) , any measurements of barriers to 
relative motion of walls of carbon nanotubes are absent up to now. Only upper limits of the shear strength 
for the relative motion of the walls along the nanotube axis were measured in a few experiments [H |2]. 
No experimental data on the barriers to relative motion of graphene layers are available. The value of the 
critical shear strength for graphite measured in the only known experiment f2T] is related to macroscopic 
structural defects of the graphite sample. Therefore the theoretical studies are particularly valuable for 
understanding phenomena and elaboration of NEMS related with relative motion and interaction of carbon 
nanotube walls. 
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For existing experimental realizations of relative motion of walls of carbon nanotubes, the movable walls 
are walls of MWNTs and have relatively great diameters: 4 — 10 nm at the relative sliding of the inner 
core of the walls attached to nanomanipulator [1], 7 — 40 nm in the nanomotors [H El [6l [7], and 3 — 7 
nm in the memory cells operating on relative sliding of the walls along the nanotube axis [8]. Therefore 
it is probably that pioneering measurements of the barriers to the relative motion of the walls will be 
performed for MWNTs with great diameter. However all previous theoretical considerations of the barriers 
to relative motion of walls deal with double-walled carbon nanotubes (DWNTs) with diameters less than 3 
nmfiaEaEallMlESESlETlIMlEHlEOlE^ That is the considered nanotubes do not correspond 

to the diameter range for which the experimental support can be expected first. Moreover the DWNTs 
with diameter greater than 5 nm have collapsed structure [33]. Therefore, calculations of the barriers 
to the relative motion of the adjacent walls of the MWNTs with a wide range of diameters, as well as 
measurements of these barriers are necessary for further progress in elaboration of nanotube-based NEMS 
and for verification of theoretical methods for consideration of interaction between nanotube walls and 
graphene layers. 

In the present work we study the dependencies of nanotube properties related with the interwall interac- 
tion on diameter of the walls. For this purpose we perform the first calculations of the interwall interaction 
energies, as well as the barriers to relative sliding of the adjacent walls along the nanotube axis for a wide 
range of the wall diameter. The results of these calculations are used to estimate the following physical 
quantities: the shear strength for the relative sliding of walls along the nanotube axis, the frequency of 
the relative axial oscillations and diffusion coefficient for the short outer wall placed in the middle of the 
nanotube, the maximum overlap of the inner and outer walls for which the controlled reversible telescoping 
can be achieved (the pulled-out core of the inner walls could be completely pushed back by capillary forces 
restoring a nanotube to its original retracted condition). The possibilities of experimental study of these 
quantities is discussed. Comparison of experimental and calculated values of these quantities can be used 
as a criterion for the correctness of the results of calculations of the barriers. 

The paper is organized as follows. Section II is devoted to analysis of the wall structure and the interwall 
interaction to prove the choice of pairs of adjacent walls for the present study. Section III presents the details 
of computational technique. Section IV gives the results of calculations of the interwall interaction energies, 
and the barriers to the relative sliding of the walls, as well as estimations of the physical quantities based 
on these barriers. Section V presents discussion of possibility of experimental study of these quantities and 
summarizes our conclusions. 

2 Structure of wall and interwall interaction 

Neglecting their structure at the ends, which can be either open or closed, carbon nanotubes are single or 
multiple layers of a cylinder rolled up from graphene sheets. Only one parameter is needed to fully determine 
the structure of the middle section of a single-walled carbon nanotube (SWNT) or a wall of a MWNT: the 
chirality index (n,m) which corresponds to a two-dimensional lattice vector c = nai -|- ma2, where ai and 
a2 are the unit vectors of graphene [35tl36j. A segment defined by the vector c becomes the circumference of 
cylindrical surface of a nanotube wall. Two types of SWNT (nonchiral SWNT) characterized by the chirality 
index of (n, n) and (ra, 0) have a simple translational symmetry, and these are referred to as armchair and 
zigzag nanotubes forming different pattern of hexagons in circumference. 

The adjacent walls of nanotube are commensurate if the ratio of the lengths of their unit cells is a 
rational fraction. In this case, a pair of walls can be considered as a quasi-one-dimensional crystal with 
the length of unit cell equal to the lowest common factor of the lengths of unit cells of constituent walls. 
If the length of the wall overlap is considerably greater than the length of the unit cell, barriers to the 
relative motion of the commensurate walls are approximately proportional to the length of the wall overlap: 
AU = AUcNc, where AUc is the barrier per unit cell and Nc is the number of unit cells corresponding to the 
wall overlap. Conceivably, it is possible to control the barriers to relative motion of the commensurate walls. 
Lack of commensurability between the adjacent walls implies a dramatic weakening of the corrugation in the 
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interwall interaction energy potential surface [22]. In contrast to the commensurate case, barriers to relative 
motion of incommensurate walls of do not increase with the length of the wall overlap, but fluctuate near 
the average value |14[ 1151 l22t I23j. These incommensurate systems, even if the wall overlap exceed micron in 
length, have barriers to the relative motion of the walls comparable to those with a wall overlap of a single 
unit cell in length. 

For adjacent commensurate walls, at least one of which is chiral, the interwall interaction energy potential 
surface is extremely flat and corrugations of the surface are smaller than the accuracy of calculations [24| 
[26| [27] . This is so because only very high harmonics of the Fourier transform of the interaction energy Ua 
between an atom of one wall and the entire other wall make contributions to the potential barriers due 
to the incompatibility of the helical symmetries of the walls [23]. Contrary to the case of chiral walls, ab 
initio calculations for the certain DWNTs ^271 |28l [29l [301 Ell Ell [33] as well as calculations using empirical 
potentials [I3l [22l [23l [Ml [SH [Ml [33] show that DWNTs with the nonchiral commensurate walls (i.e. both 
armchair or both zigzag) have considerable barriers to relative motion of walls. Therefore, the pairs of 
the nonchiral commensurate walls can be considered as potential components in NEMS for which a precise 
control of motion or position of the walls is required. A wide set of such NEMS which operating principle 
are based on the dependence of the conductivity or /and the interwall interaction energy on the relative 
position of the walls on subnanometer scales was proposed. This set includes now a varying nanoresistor 
[141 [T5l [37l [38] , a stress nanosensor [39] , an electromechanical nanothermometer [401 |4T] , a nanoactuator for 
the transformation of a force directed along the nanotube axis to the relative rotation of its walls [4T] 142] . 
and a nanoresonator |43j . Thus, analysis of the relative motion and interaction of nonchiral commensurate 
carbon nanotube walls holds the key to success of these applications. Moreover, due to considerable value of 
the barriers to relative motion of walls, nanotubes with such walls should be considered as best candidates 
for measurements of these barriers. 

3 Computational technique 

Since the system consisting of two nonchiral commensurate walls with great radius has the unit cell which 
contains up to several thousands atoms, ab initio calculations are not available for such a system. Thus we 
perform empirical calculations of the interwall interaction energy and the barriers to the relative motion of 
the nonchiral commensurate walls with the aim to obtain the qualitative dependencies of these quantities 
on the wall radius. The interaction between atoms of adjacent walls is described here by 6-12 Lenard-Jones 
potential U = 4e(((T/r)^^ — (a/r)^). Note that contribution of many-body terms into the complete Van 
der Waals interaction is not essential at relative motion of linear clusters interacting along its full length 
[44]. The parameters of potential, e = 2.968 meV and a = 3.407 A, were fitted to interlayer distance and 
modulus of elasticity of graphite and used in study of the interwall interaction energypotential surface of the 
DWNTs [m [Ml [HI [26] . The upper cutoff distance of the potential is chosen so that calculate the barriers 
to the relative sliding of the walls with the accuracy 0.1 %. Namely, the upper cutoff distance is 42, 120, 25 
and 35 A for (n,n)@(n+5,n+5), (n,n)@(n+6,n+6), (n,0)@(n+9,0), and (n,0)@(n+10,0) pairs of adjacent 
walls, respectively. The length of the outer wall is equal to the length of the unit cell and the length of 
the inner wall is chosen so that all pairs of atoms with interatomic distances within the cutoff distance are 
taken into consideration. The structure of the walls is obtained by mapping of graphene with bond length 
Co = 1.41 A on cylindrical surface. The walls are considered as rigid. Account of the deformation of walls is 
not essential for the shape of potential surface both for the interwall interaction energy of DWNTs [22] and 
the intershell interaction energy of carbon nanoparticles [45 i 146] . For example, the barriers to the relative 
wall rotation and sliding along the axis for the (5,5)@(10,10) DWNT calculated for walls with unannealed 
structure |26] coincide within 14 % with results of Dresselhaus et al. (used annealed wall structure) |13] . 
Note also that 6-12 Lenard-Jones potential was used to study the dependencies of the ground state interwall 
interaction energy [471 SH] and radial breathing mode frequencies [48] as function of radius and chiral angle 
of walls for DWNTs with chiral walls and to determine characteristics of nanotube-based bolt/nut pairs |49] . 
According to the classification scheme of DWNTs with commensurate wall [26] (which is also the clas- 
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Figure 1: The schematic diagram in cyhndrical coordinates for the relative positions of walls of the 
(5,5)@(10,10) DWNT corresponding the maximum of the interwall interaction energy. Bonds and atoms of 
the inner (5,5) wall are shown by thin lines and grey circles, respectively. Bonds and atoms of the outer 
(10,10) wall are shown by thick lines and black circles, respectively. Point A corresponds to the coincident 
symmetry axes U2 of the inner and outer walls which passes through the midpoint of the carbon bond, point 
B corresponds to the coincident axis U2 of the inner wall which passes through the midpoint of the carbon 
bond and axis U2 of the outer wall which passes through the center of the hexagon, point C corresponds to 
the coincident axis U2 of the inner wall which passes through the center of the hexagon and axis U2 of the 
outer wall which passes through the midpoint of the carbon bond, point D corresponds to the coincident axes 
U2 of the inner and outer walls which passes through the center of hexagon. Vector a = {6z/2,6^/2) shows 
the displacement and rotation of one of the walls to the relative position corresponding to the minimum of 
the interwall interaction energy. 
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sification scheme of pairs of coaxial commensurate walls) such pairs form families with the same interwall 
distance and unit cell length for each pair of any family. Here we consider here four families of such pairs: 
(n,n)@(n+5,n+5), (n,n)@(n+6,n+6), (n,0)@(n+9,0), and (n,0)@(n+10,0) with the interwah distances 3.37, 
4.04, 3.50 and 3.89 A, respectively. The inner wall chirality index n up to 250 are considered, such chirality 
index corresponds to the outer wall radiuses from 0.65 to 17 nm for the pairs of armchair walls and from 
0.7 to 10 nm for the pairs of zigzag walls. 

To analyze the nanomechanical properties connected with the relative motion and interaction of nanotube 
walls, we must calculate the ntershell interaction energy potential surface, i.e. the dependence of the energy 
U of the interaction between two adjacent walls on the coordinates describing the relative position of the 
walls (angle (p of the relative rotation of the walls about the nanotube axis and length z of the relative 
displacement of the walls along this axis). 

The relative positions of the walls corresponding to a higher symmetry of the system are critical points 
of interwall interaction energy U{z,(j)) (extrema or saddle points) [50]. For such relative positions of the 
walls, some of the second-order axes U2 of the inner and outer walls are in line. The second-order U2 
axis is perpendicular to the principal axis of the wall and passes through the midpoint of the carbon bond 
or the center of the hexagon. The detailed consideration for the case of coaxial nonchiral commensurate 
walls is presented in [25\ I31j . In particular, the possible relative positions corresponding to the critical 
points are listed in [31], and the relative positions corresponding to the minima of the interwall interaction 
energy are determined using 6-12 Lenard- Jones potential in ^2^. In accordance with the symmetry of 
armchair (ni, ni)@(n2, ^2) and zigzag (ni, 0)@(n2, 0) DWNTs, the unit cell of the interwall interaction 
energy potential surface of such pairs of walls is a rectangle with sides 5z = ta/2 and 6^ = iTN/nin2, where 
ta is the length of the unit cell of the wall and is the greatest common division of numbers ni and n2 
[50) . As an example we present in figure 1 the schematic diagram describing the relative positions of walls 
of the (5,5)@(10,10) DWNT corresponding to the maximum of the interwall interaction energy. 

Empirical calculations with the use of the Lennard-Jones potential show that for all considered armchair 
{n,n)@{n+5,n+5) DWNT with n = 5 - 15 and zigzag (n,0)@(n-h9,0) DWNT with n = 9 - 18, the interwaU 
interaction energy can be successfully interpolated within an accuracy of about 1% of the value of the 
barriers by the following expression [26] : 




where Uq is the mean energy of interwall interaction, AU^ and AC/^ are the barriers to the relative sliding 
of the walls along the nanotube axis and the relative rotation of the walls (between their relative positions 
corresponding to the interaction energy minimum), respectively. Henceforth, we will normalize the values of 
Uq, and AUz to an atom of the outer wall. This empirical results is confirmed for several DWNTs with the 
use of density functional method within the accuracy of the calculations \27 \ \32\ \5T\. In this case, where the 
expression ([T|) is adequate, the four types of the critical points on the surface are the global minimum, global 
maximum, and two saddle points. Thus to determine the barriers to the relative motion of the walls the 
interwall interaction energy is calculated here at the critical points. For DWNTs with compatible rotational 
symmetry of the walls (A^ 7^ 1), the barrier to the relative rotation, AU^, can be significant. For DWNTs 
with incompatible rotational symmetry of the walls {N = 1), both ab initio [31] and empirical [24:\ [26] 
results show that the dependence of the interwall interaction energy on the angle (p is extremely flat, thus 
the second term in expansion ([T]) can be neglected. The present work is devoted to pairs of adjacent walls 
with great radius. The rotational symmetry of the walls is incompatible for all such pairs. Therefore only 
the barriers to the relative sliding of the walls along the nanotube axis and physical quantities related with 
these barriers are calculated. 
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Figure 2: Calculated interwall interaction energy (in meV per atom of the outer wall) as a function of the 
radius of the outer wall (in nm) of the (n,n)@(n+5,n+5) (red solid thick line), (n,n)@(n+6,n+6) (blue 
dotted thing line), (n,0)@(n+9,0) (green dashed thick line), and (n,0)@(n+10,0) (black solid thing line) 
families of wall pairs. 

4 Properties related with interwall interaction 

The calculated dependencies of the interwall interaction energy corresponding to the minimum of the function 
U{z,(j)) on the radius of the outer wall are presented in figure 2 for four mentioned above families of the 
adjacent wall pairs. These values of the interwall interaction energy lie in the range 24 — 50 meV per 
atom of the outer wall. This result is in good agreement with the experimental values of the energy of 
interaction between graphite layers 52ib5 meV/atom obtained recently from the experiments on the thermal 
desorption of polyaromatic hydrocarbons from a graphitic surface fTSi , 43 meV/atom from the heat-of- 
wetting experiment |20] . 35ibl0 meV/atom from an analysis of the structure of collapsed carbon nanotubes 
[18], and between the walls of MWNTs 20—33 meV/atom [2]. The figure 2 shows that for pairs from the 
same family the interwall interaction energy tends to a constant value with increase of the outer wall radius. 
We found that for the great radiuses of the outer walls the dependencies of the interwall interaction energy 
on the outer wall radius can be interpolated by the following expression: 
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C/ = C/oo + f/iexp(^-^) , (2) 

where Uoo, Ui and R' are the fitting parameters, R2 is the radius of the outer wall. The calculated values 
of the fitting parameters are listed in table 1. Note that the fitting parameter R' , which characterizes the 
rate of exponential decay of the interpolated dependencies, coincides within the calculation accuracy for 
both considered families of the armchair adjacent wall pairs of and for both considered families of the zigzag 
adjacent wall pairs. The parameter Ui is order of magnitude less than the parameter Uoo for all considered 
families of the pairs of the adjacent wall. Therefore the interwall interaction energy per one atom of the 
outer wall only slightly increases with the radius increase and tends to a constant value Uoo for these pairs 
of walls with the radius greater than 5 nm. 

Table 1 Calculated fitting parameters Uoo, Ui and R' of the interpolation by expression ([2D of the 
dependencies of the interwall interaction energy on the radius of the outer wall. The interpolation is 
performed for the radius greater than 5 nm. 



family of wall pairs 


Uoo (meV/atom) 


Ui (meV/atom) 


R' (nm) 


(n,n)@(n+5,n+5) 
(n,n)@(n+6,n+6) 
(n,0)@(n+9,0) 
(n,0)@(n+10,0) 


-49.713±0.004 
-35.944±0.003 
-49.132±0.003 
-39.942±0.004 


3.65±0.03 
3.18±0.03 
4.87±0.03 
4.43±0.02 


4.61±0.04 
4.64±0.04 
3.40±0.02 
3.39±0.02 



We confirm also the previous qualitative results of atomistic empirical [13\ [26l W7\ HH] and ab initio 
calculations [31] for DWNTs with the smaller radiuses: 1) the interwall interaction energy per one atom 
increases with an increase of the radius of DWNT [26[ I31j : 2) the interwall interaction energy per one 
atom is greater for DWNTs with the interwall distance close to the distance between graphite layers 3.4 

A [la [311 in SHI- 

The analogous tendency to a constant of the interwall interaction energy per one atom with the increase 
of wall radius was obtained in the framework of approach with continuous distribution of atom density not 
only for interaction between pairs of adjacent walls but also for the case where the interaction with the next 
walls of the inner core or outer shell of the MWNT is taken into account [Tl]. That is the account of this 
additional interaction energy leads only to increase of the parameter Uoo and this increase is found to be 
within 27 % for any number of walls in the MWNT pjj. Such a change of the parameter Uoo is comparable 
to the scatter in the listed above experimental data for graphite interlayer interaction energy. Therefore we 
believe that the corrections to the interaction energy due to the interaction with the next walls of the inner 
core or outer shell of the MWNT obtained in the framework of the approach with continuous distribution 
of atom density pTj can be applicable also to the our the results obtained in the framework of the atomistic 
approach within the accuracy of both approaches. Thus the atomistic calculations of the mentioned above 
corrections to the interaction energy is beyond the purpose of the present work which is devoted mainly to 
the physical properties related with the barriers to the wall motion. 

A telescopic extension of the inner wall gives rise of a capillary force Fc pulling the inner wall back into 
the outer wall. The average value of this force is defined as 

dU \ _ Uo4n2 , . 

rU / ~ f ' ^ ' 

where Lo„ is the overlap length of the walls, An^ is the number of atoms in the unit cell of the outer wall. 
The calculated dependencies of the average capillary force on the radius of the outer wall are presented in 
figure 3. Since the interwall interaction energy tends to a constant value for the pairs of walls from the 
same family with the increase of the wall radius, the average capillary force for such pairs is approximately 
proportional to the radius. 
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Figure 3: Calculated capillary force pulling the inner wall back into the outer wall (in nN) as a function of 
the radius of the outer wall (in nm) of the (n,n)@(n+5,n+5) (red solid thick line), (n,n)@(n+6,n+6) (blue 
dotted thing line), (n,0)@(n+9,0) (green dashed thick line), and (n,0)@(n+10,0) (black solid thing line) 
families of wall pairs. 

The calculated dependencies of the barriers AC/^ for the relative sliding of walls on the outer wall 
radius are shown in figure 4 for all considered families of the wall pairs. The relatively small increase of the 
interwall distance leads to decrease of of the barrier to the relative sliding of the walls by orders of magnitude 
(compare (n,n)@(n+5,n+5) and (n,n)@(n+6,n+6), (n,0)@(n+9,0) and (n,0)@(n+10,0) families). Therefore 
the contribution of the interaction between the other walls of the MWNT into the barrier to the relative 
sliding of the walls is negligible in comparison with the contribution of the interaction between the adjacent 
walls. Note that the barrier to relative sliding of walls of (5,5)@(10,10) pair of walls calculated here is close 
to that calculated using density functional method [32]. Analogously to the dependencies of the interwall 
interaction energy on the outer wall radius, the dependencies of the barriers to relative sliding of walls for 
the great radius of the outer wall are interpolated by the following expression 



where AUoo, AUi and i?^ are the fitting parameters. The calculated values of the fitting parameters for 
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Figure 4: Calculated barrier to the relative sliding of the walls (in meV per atom of the outer wall) as 
a function of the radius of the outer wall (in nm). (a) the (n,n)@(n+6,n+6) family of wall pairs (green 
triangles) and the (n,n)@(n+5,n+5) family of wall pairs (red circles), (b) the (n,0)@(n+10,0) family of wall 
pairs (green triangles) and the (n,0)@(n+9,0) family of wall pairs (red circles). 



families of wall pairs with greater barriers (n,n)@(n+5,n+5), (n,0)@(n+9,0), and (n,0)@(n+10,0) are listed 
in table 2. Analogously to the dependencies of the interwall interaction energy on the radius of the outer 
wall, the parameter AC/i is order of magnitude less than the parameter AUoo for all families of the wall pairs 
presented in table 2. Therefore for these pairs with the outer wall radius greater 5 nm interwall interaction 
energy per one atom of the outer wall only slightly increases with the radius increase and tends to a constant 
value AC/qo- 



Table 2 Calculated fitting parameters AUoo, AC/i and i?^ of the interpolation by expression of the 
dependencies of the barrier to relative sliding of the walls on the radius of the outer wall. The interpolation 
is performed for the radius greater than 5 nm. 



family of wall pairs 


AC/oo (meV/atom) 


AUi (meV/atom) 


i?^ (nm) 


(n,n)@(n+5,n+5) 

(n,0)@(n+9,0) 

(n,0)@(n+10,0) 


0.0104 
0.237 
0.0330 


0.0010 
0.013 
0.0020 


3.48±0.18 
5.83±0.03 
5.95±0.07 



The calculated values of barriers AUz for the relative sliding of the adjacent walls can be used for 
calculating the values of a number of physical quantities characterizing the interaction and relative motion 
of the nanotube walls. Expression ([T|) for the dependence of the interwall interaction energy U{(p, z) on the 
relative position of the walls defines the threshold static friction force Fz for the relative sliding of the walls 
along the nanotube axis: 



An2'KLovAUz 

Ozta 

The shear strength for relative sliding of the walls along the nanotube axis is defined as 
where S is the area of the overlap surface of the walls. 



(5) 



(6) 
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S = 27rLo. [ ^^^^^ )= TTLo^iRi + R2), (7) 



2 

where i?i is the radius of the inner wall. Figure 5 presents the dependencies of the shear strength on the 
radius of the outer wall for four considered families of the wall pairs. 
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Figure 5: Calculated shear strength for relative sliding of the walls (in MPa) as a function of the radius of the 
outer wall (in nm). (a) the (n,n)@(n+6,n+6) family of wall pairs (green triangles) and the (n,n)@(n+5,n+5) 
family of wall pairs (red circles), (b) the (n,0)@(n+10,0) family of wall pairs (green triangles) and the 
(n,0)@(n+9,0) family of wall pairs (red circles). 

The equality of forces Fc = acting on the movable wall controls the minimal overlap length Lm of the 
walls for which the static friction force F^ prevents the inner wall from being pulled by the capillary force 
Fc- Substituting Eqs. ([3]) and ([5]) for forces Fc and Fz, respectively, in this equality we have 

= (8) 

For the (n,n)@(n+6,n+6) family the calculated value of the minimal overlap length is about 14 /xm 
which is too much for nanotubes produced. The dependencies of the minimal overlap length Lm on the 
radius R2 of the outer wall for other considered families of the wall pairs are presented in figure 6. 

Let us consider the case where the outer wall is shorter than adjacent wall and placed near middle of 
the adjacent wall without telescopic extension so that the distance between the edges of the outer wall and 
the adjacent wall is considerably greater than the parameter a of the potential. In such a case influence of 
the wall edges on the interwall interaction energy can be disregarded. To find the frequency of small axial 
oscillations of the outer wall in this case, we approximate expression ([T|) for the dependence of the interwall 
interaction energy on their relative position by a parabolic potential well in the vicinity of the minimum Um 
of this dependence, 

U{z') = Um + ^, (9) 

where z' is the wall displacement along the nanotube axis relative to the energy minimum of interwall 
interaction. For (ni, ni)@(n2, 71-2) and (ni, 0)@(n2, 0) pairs of walls, we have 

fe = ?=H^(|)'. (10) 
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Figure 6: Calculated minimal overlap length of the walls for which the static friction force prevents the core 
of the inner walls from being pulled by the capillary force (in nm) as a function of the radius of the outer 
wall (in nm). (a) the (n,n)@(n+5,n+5) family of wall pairs (circles), (b) the (n,0)@(n+10,0) family of wall 
pairs (green triangles) and the (n,0)@(n+9,0) family of wall pairs (red circles). 



Frequency Vz of small relative axial oscillations of the wall is defined as 




(11) 

where M is the mass of the outer wall, M = An2mQL/ta, where L is the length of the outer wall. The 
expression for the frequency of the axial oscillations of the outer wall assumes the form 



'2-^. (12) 

mo 

Figure 7 shows the dependencies of this frequency on the radius of the outer wall for four considered families 
of wall pairs. 

If the energy of thermal motion of a short outer wall is much lower than the barrier At/ to the sliding of 
this wall along the nanotube axis [kT <C AC/, where AC/ = 4:n2AUzL/ta), then hopping diffusion or drift of 
the shorter wall along the nanotube axis takes place. Recently such a drift actuated by imposing of thermal 
gradient was observed for short shell consisting of several outer walls of many- walled carbon nanotube [6]. 
The Fokker-Planck equation has been obtained in a general case of diffusion and drift of the walls along the 
helical line |14^ [T5] . In the case of the motion of the mobile wall along the nanotube axis diffusion coefficient 
Dz is defined as 

Ir^ ,2 _____ f 4n2AC/^L 



°' = 5"'«-H-^f^rJ (13) 

where 0^ is the preexponential factor in the Arrhenius formula for the frequency of the jump of the mobile 
wall between equivalent minima of the interwall interaction energy. 

Simulation based on the molecular dynamics method shows that the preexponential factor in the Ar- 
rhenius formula for the frequency of the jump of a C6o@C24o nanoparticle shell between equivalent minima 
of intershell interaction energy is = 540 it 180 GHz [l5]. We found that frequency of small relative rota- 
tional oscillations near the minima of intershell interaction energy calculated here for such a nanoparticle 
is = 60 GHz for the same shape of the shells and the same potential of interaction between atoms of 
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Figure 7: Calculated frequencies of relative axial oscillations of the outer wall (in GHz) as a function of 
the radius of the outer wall (in nm). (a) the (n,n)@(n+5,n+5) family of wall pairs (green circles), the 
(n,n)@(n+6,n+6) family of wall pairs (red squares); (b) the (n,0)@(n+9,0) family of wall pairs (green 
circles), the (n,0)@(re+10,0) family of wall pairs (red squares). 



the shells. Thus, ratio Vt/v the preexponential factor in the Arrhenius formula to the frequency of small 
oscillations is approximately 10. We believe that ratio ^/v for other carbon nanostructures with enclosed 
graphene layers (in particular, carbon nanotubes) is on the same order of magnitude. We used the value 
of ^zl'^z ~ 10 for estimating the diffusion coefficients for diffusion of the outer wall along the nanotube 
axis. For the (n,n)@(n+6,n+6) family of the wall pairs the ratio ISJJ jkT > 0.1 for the short outer wall of 
length 100 nm for all considered pairs of walls. Thus the expression (113p is not adequate for estimation of 
the diffusion coefficient D^. for this family of the wall pairs. The calculated dependencies of the diffusion 
coefficients for the motion of the short outer wall 10 nm in length moving along the infinite inner wall on 
the radius of the outer wall are presented in figure 8 for other three considered families of the wall pairs. 
Only pairs with AU/kT < 0.1 are shown. The lower limit Dz = 10~^^ nm^/ns of the diffusion coefficient 
range shown in figure 8 corresponds to displacement of the outer wall by ~ 1 nm during a day. 

Under the action of force directed along the nanotube axis, the short outer wall may drift along the 
axis at a velocity v = BzFd, where Bz is the mobility of the outer wall for its motion along the axis. In the 
case of a small force, F^Sz <C 2AU, mobility Bz is related to the diffusion coefficient by the Einstein relation 
Dz = kTBz [Hill]. 

5 Discussion and conclusions 

The dependencies of interwall interaction energies, as well as barriers to relative sliding of the walls along 
the axis, on the radius of the outer wall are calculated for the pairs of adjacent nonchiral commensurate 
walls of carbon nanotubes with the wide range of radiuses. It is found that for walls with radius greater 
than 5 nm both the interwall interaction energy and barriers to relative sliding of the walls per one atom of 
the outer wall only slightly increase with the radius increase and tend to the constant values. The results 
of these calculations are used to obtain for these pairs of walls the dependencies on the radius of the outer 
wall for the wide set of the measurable physical quantities related to the calculated barriers. Certain of 
the calculated physical quantities which are determined only by the barrier to sliding per one atom of the 
outer wall and do not depend on the radius of the outer wall also only slightly increase with the radius 
increase and tend to the constant values. Such quantities are the shear strength for the relative sliding of 
the walls along the nanotube axis, the frequency of the axial oscillations of the outer wall and the maximum 
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Figure 8: Calculated diffusion coefficient of outer wall of length 10 nm (in nm^/ns) as a function of the 
radius of the outer wall (in nm) of the (n,n)@(n+5,n+5) (red solid thick line), n,0)@(n+9,0) (green dashed 
thick line), and (n,0)@(n+10,0) (black solid thing line) families of wall pairs. 

overlap of the inner core and outer shell of the walls for which the controlled reversible telescoping can be 
achieved (the pulled-out core could be completely pushed back by capillary forces restoring the nanotube to 
its original retracted condition). The diffusion coefficients for the relative sliding of the walls along the axis 
of the nanotube is determined by the total barrier to the relative sliding and exponentially decrease with 
the radius increase. 

Let us consider how to experimentally verify our results. At present, only the upper boundary of the shear 
strength for the relative sliding of the walls along the nanotube axis is determined with the help of atomic 
force microscopy is < 0.04 MPa [2]. In most cases, the barriers to the sliding of adjacent walls (and, 
hence, the corresponding shear strength) are negligibly small (due to the incompatibility of translational 
symmetry of the walls for incommensurate walls [23] and to incompatibility of helical symmetries of the walls 
for commensurate chiral walls [22l [26| [27] ) . In experiment [2] , the shear strength was measured only for a 
single pair of walls of a MWNT, which was characterized by a smaller value of this strength as compared 
to other pairs of adjacent walls. The chirality indices of the walls were not determined. Thus, we believe 
that the experimental value of the upper boundary of the shear strength for the relative sliding of the walls 
corresponds to incommensurate or commensurate chiral walls. However, the barrier to the relative sliding of 
adjacent walls may attain an appreciable value only for commensurate nonchiral walls [2 ^ \23 \ IM ) 126). The 
shear strength calculated here for such pairs of walls exceed by several orders of magnitude the corresponding 
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boundary of this value (0.04 MPa) measured by an atomic force microscope for the relative sliding of adjacent 
walls with indeterminate chirality indices. Consequently, the shear strength for these pairs of walls can easily 
be determined experimentally. In our opinion, the minimal overlap length of the walls in the case of the 
telescopic extension of the inner core of the walls, for which the static friction force prevents the inner core 
from being pulled by the capillary force, can also be measured using atomic force microscopy. 

The nanotube with a and short outer wall can be obtained via the electrical-breakdown technique j52tl53j. 
which consists of passing a large current through the nanotube and was used for fabrication of nanotube- 
based nanomotors O EJ [7]. Such a DWNT with the short outer wall can be used for measurements of the 
frequency of axial oscillation of the outer wall relative the fixed inner core of wall using terahertz spectroscopy 
or Raman spectroscopy. It should be noted that Raman spectra have recently been obtained for isolated 
nanotubes |54J . The diffusion coefficient for the short outer wall moving along the nanotube axis can be 
measured with the help of scanning electron microscopy which was used to study the motion of the short 
outer wall driven by thermal gradient [6]. 

Thus we believe that the calculated physical quantities which are determined by the barriers to relative 
sliding of walls along nanotube axis can be measured by contemporary experimental techniques. Such 
measurements will give the first results for verification of the theoretical methods for studying the interwall 
interaction in nanotubes and would facilitate progress in developing nanoelectromechanical systems based 
on the relative motion of the walls of nanotubes. 
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